Dynamics of Multiferroic Domain Wall in Spin-Cycloidal Ferroelectric DyMnOs 
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We report the dielectric dispersion of the giant magnetocapacitance (GMC) in multiferroic 
DyMnOa over a wide frequency range. The GMC is found to be attributable not to the soft- 
ened electromagnon but to the electric- field-driven motion of multiferroic domain wall (DW). In 
contrast to conventional ferroelectric DWs, the present multiferroic DW motion holds extremely 
high relaxation rate of ~10^ s~^ even at low temperatures. This mobile nature as well as the model 
simulation suggests that the multiferroic DW is not atomically thin as in ferroelectrics but thick, 
reflecting its magnetic origin. 

PACS numbers: 75.47.Lx, 75.80.+q, 75.60.Ch 



Prospective materials toward electric control of mag- 
netism in a solid are multiferroics, in which electric and 
magnetic orders coexist P, 0] ■ In conventional ferromag- 
netic or ferroelectric materials, the motion of ferroic do- 
main walls (DWs) is a key to the functions; it provides 
huge linear-response (i.e., magnetic or dielectric suscep- 
tibility) as well as low-field control of the ferroic 
domain. By analogy, the control of the multiferroic [i.e. 
concurrently (anti)ferromagnetic and ferroelectric] DWs 
in multiferroics may provide a new prospective means 
to attain the electric (magnetic) control of the magnetic 
(ferroelectric) domain on a macroscopic scale as well as 
the enhancement of the dynamical magnetoelectric (ME) 
susceptibility. In this Letter, through dielectric measure- 
ments over a wide frequency region on the multiferroic 
perovskite, DyMnOs [y], we identify the microscopic mo- 
tion of multiferroic DW as a source of large magneto- 
electric coupling. The multiferroic DW is found to be 
dynamically active even at low temperatures, enabling 
the electric control of magnetic domains via the macro- 
scopic DW movement. These characteristics suggest that 
the present multiferroic DW is thick in width, in contrast 
to the conventional thin ferroelectric DW. 

In multiferroic materials, ferroelectric polarization P 
is often induced by magnetic order through the inverse 
Dzyaloshinskii-Moriya (DM) mechanism [7|. In the per- 
ovskite manganites /JMnOa {R = Tb and Dy) with 
GdFeOs-type distortion [Fig. 1 (a)] , there is growing ev- 
idence that the cycloidal magnetic order induces the fer- 
roelectricity through the inverse DM mechanism [1, 
i.e., P oc J2i j ^ij ^ i'^i ^ "^i)' where Si {Sj) is the 
electron spin vector at site i (j) and Cjj is the unit vector 
connecting the two sites. In TbMnOs or DyMnOa, the 
&c-cycloidal spin order is emergent at low temperatures 
with the propagation vector q along the b axis [q\\b, see 
Fig. 1(b)] [1,1a 113]; thus P is induced along the c axis. 
It is also confirmed that a magnetic field H\\b induces 



the polarization flop from -Pjjc to ■ Pjja with the direction 
of q (II 6) unchanged [H, [III, [H, [Tj], as exemphfied in 
Figs. 1(d) and 1(e) for DyMnOa. This behavior is at- 
tributable to the flop of the spin cycloid plane from be to 
ab [Fig. 1(c)]. The salient ME feature characteristic of 
DyMnOs, which is investigated here, is the giant magne- 
tocapacitance (GMC) emerging concomitantly with the 
magnetic- field- induced P fiop [6|: as shown in Fig. 1(f), 
the dielectric constant (-E||a) shows a large enhancement, 
from ~40 to ~130 (in the case of 10 K), in the course of 
the magnetic-field-induced transition. This phenomenon 
may give an important clue to the general strategy to 
realize the colossal ME effect. 

We have investigated the magnetocapacitance with 
varying frequencies of electric field (-B||a) from 1 KHz to 
500 MHz using two equipments: a LCR meter (Agilent 
E4980A) for 1 KHz-2 MHz and an impedance analyzer 
(Agilent E4991A) for 1 MHz-500 MHz. Figure 1(f) dis- 
plays the magnetic-field dependence of the real part of di- 
electric constant, ea' (£^||a). Note that the GMC emerges 
below 100 KHz but is considerably suppressed above 10 
MHz. This behavior clearly demonstrates that the GMC 
has prominent dielectric dispersion around 1 MHz at 10 
K. More insight into the origin of GMC is obtained from 
the spectral analysis. We measured the spectra of real 
and imaginary parts of dielectric constant, Ca' and Ca" 
(£^||a), at various magnetic fields and temperatures. As 
a typical example, the results at 10 K are shown in Figs. 
2(a) and 2(b). These spectra ensure that the GMC is 
a phenomenon emergent only below 10^-10^ Hz. A new 
aspect revealed here is that the spectral shape of GMC 
is not the resonance type but the relaxation type, indi- 
cating that the GMC does not arise from bosonic exci- 
tations, such as soft-mode phonon s 1141 a nd electrically 
active magnons (electromagnons) |15l . Ilq . This finding 
is clearly incompatible with the prevailing argument that 
the softening of electromagnon results in the GMC at low 
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FIG. 1; (Color online) (a) Crystal structure of DyMnOs. 
(b) and (c) Magnetic structures of &c-cycloidal spin order 
(b) and afo-cycloidal spin order (c). (d) Magnetic field- 
temperature phase diagram of DyMnOg for the case of H\h. 
The phase boundaries were determined from the dielectric 
measurements. The closed squares and circles represent the 
transition points in increasing and decreasing field processes, 
respectively, (e) Magnetic-field dependence (i?||6) of polar- 
ization at 10 K along the c axis (P||c) and the a axis (P||a). 
(f) Magnetic-field dependence (i?||6) of dielectric constant at 
10 K along the a axis at various frequencies. 



frequencies [l^, The origin of dielectric relaxation 
can be deduced by considering several features: (i) the 
magnetic-field dependence of the magnitude of relaxation 
mode, Aea', (ii) the temperature dependence of the relax- 
ation rate at the flop transition, and (iii) the value of the 
relaxation rate itself. Here we evaluated Aea'(i7, T) by 
ea'(i?, T, IKHz) - ea'(i?, T, 500MHz) and defined the re- 
laxation rate as 1/t{H,T) = 2T:fpca.k{H,T), where /peak 
is the frequency of Eq" peak. Below we refer to (i)-(iii) 
and conclude that the GMC is provided by the motion of 
multiferroic DWs between the f ||±a (&c-cycloidal) and 
the f ||±c (a6-cycloidal) domains. Hereafter, the relevant 
multiferroic DWs are abbreviated as 'DW±a/±c [see Fig. 
3(a)]. 

The magnetic-field dependence of Aeq' and l/r at var- 
ious temperatures are shown in Figs. 2(c) and 2(d), re- 
spectively. Note that Acq' shows a sharp peak even at 
5 K. This tendency is not expected for the case that the 
GMC were due to the polarization fluctuations; at low 
temperatures the polarization flop transition is of the first 
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FIG. 2: (Color online) (a) and (b) Spectra of dielectric con- 
stant at 10 K under various magnetic fields: real part (a) and 
imaginary part (b). (c) and (d) Magnetic-field dependence of 
the magnitude of the relaxation mode (c) and the relaxation 
rate (d) at given temperatures. Inset in (d): Temperature 
dependence of the relaxation rate at the fiop transition, I/tc 
(double logarithmic plot). The solid line is a guide to the eye. 



order and thus the enhancement of fluctuations is un- 
likely. In Fig. 2(d), we focus on the local minimum value 
of \/t[H) at a given temperature, say 1/tc (i.e., l/r at 
the polarization flop transition) . Note that although the 
1/t relevant to polarization fluctuations should become 
larger toward low temperatures in general the 
currently observed 1/tc rather diminishes. Therefore, in 
terms of (i) and (ii), it is unlikely that the observed re- 
laxation comes from the polarization fluctuations. As the 
origin of such a relaxation characteristic, the microscopic 
displacement of DW under oscillating E is the most plau- 
sible. 

Among the various possible multiferroic DWs in 
DyMnOa [Fig. 3(a)], the relevant ones that may give 
rise to large ej response are the motion of lyW ^a/~a 
and/or ±a/±c- In the P||a phase, we can identify the 
DW-|.a/-a contribution to from the difference between 
Eq' in i-'||±a multi-domain state and in P||+a single- 
domain state. We found that although the DW_|_a/_a 
contribution is quite small (Aea' < 7), the WNj^a/-a i™- 
tion also shows the relaxation-type spectra and its l/r 
is comparable to 1/rc of GMC (not shown). Therefore, 
from the quantitative point of view, the GMC emerging 
only near the flop transition is also attributable to one of 
the DW motions, namely the iyW±a/±c motion [l^. Re- 
markably, this DW±a/±c motion holds a high relaxation 
rate of ~ 10'' s~^ even at low temperatures, in contrast to 
the conventional ferroelectric DW motion, which usually 
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FIG. 3: (Color online) (a) Various multiferroic DWs conceiv- 
able in DyMnOs. The red lines represent the DW through 
which both the spin helicity {Si x Sj) and the polariza- 
tion (P) are reversed, while the blue lines represent the DW 
through which they rotate by 90 degrees, (b) Initial P-E 
curves under various magnetic fields at 5 K, demon- 

strating electric-field-induced movement of multiferroic DW 
between the P||c (6c-cycloidal) and the P||a (ab-cycloidal) 
domains. 



freezes at low temperatures [20|. In fact, the tempera- 
ture dependence of I/tc obeys a power-law rather than 
the Arrhenius law [the inset of Fig. 2(d)], resulting in the 
gradual slowing down of 1 /tc toward low temperatures 
and thus in the high relaxation rate even at 5 K. 



Since the DW. 



±a/±. 



motion is still unfrozen even at the 



lowest temperature, 5 K, the macroscopic movement of 
depinned DW-|-a/-|-c is expected under a strong electric 
field. In conventional ferroelectrics, if one attempted to 
induce the macroscopic movement of DW at 5 K, an im- 
practically strong electric field or long waiting time would 
be required because the DW motion via nucleation on the 
DW is hard to occur at such low temperatures 
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In DyMnOa, however, the macroscopic DW±a/±c move- 
ment does occur even at 5 K; figure 3(b) displays the 
"initial" P-E curve (-E||a) under various magnetic fields. 
Note that the hysteresis behavior is observed only near 
the flop transition field (~ 1.71 T at 5 K). This hystere- 
sis indicates that the initial and final states are different. 
Since we took the experimental procedure to suppress 
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FIG. 4: (Color online) Calculated DW structure between the 
P||-|-c (6c-cycloidal) and P||+a (a6-cycloidal) domains for 36 
X 6 Mn sites. Blue and red arrows represent the Mn spins 
and the local electric polarizations, respectively. The color 
gradation represents the angle of local electric polarization 
relative to the a axis. The angle becomes 45 degrees along 
the DW center, which runs parallel to the 6 axis. 



the mixture of P\\-a domain in the P|la phase [23|, this 
irreversible increase in Pa can be regarded as the conse- 
quence of the macroscopic movement of Vi^ j^aj+c and/or 
DW_|_(j/_c- Therefore, the irreversible P-E curve suggests 
that the magnetic cycloidal plane was changed from he to 
ah by an electric stimulation through the DW movement. 



The fact that the DW. 



±a/±c 



motion does not freeze 



even at 5 K implies that the DW-to/j-c motion does not 
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occur via the nucleation process on the DW 
Moreover, it is known that there is a close relationship 



between the DW mobility and the DW thickness [24 1. 



In this context, the DW-ta/j-c structure is of particu- 
lar interest and importance. Here we postulate that the 
DW structure is determined by minimizing mainly the 
magnetic energy cost, because the ferroelectricity of mag- 
netic origin in DyMnOa is quite weak in magnitude com- 
pared with conventional ferroelectrics. To calculate the 
internal structure of DW+a/+c, we employed the two- 
dimensional classical J\-Ji model with ferromagnetic J\ 
(= 0.8 meV) on the tetragonal x and y bonds and antifer- 
romagnetic Ji (= 0.96 meV [2^) on the diagonal bonds 
along the h axis: these values reproduce a spiral spin or- 
der with q = 0.36 (||6) @. To mimic the experimental 
situation where the ah- and 6c-plane spin cycloids are de- 
generate at the phase boundary, we added the single-ion 
anisotropy term, DY,^{Sl^ + S^^) with D = 0.20 meV, 
and the orthorhombic anisotropy term, f3J2i^ai^ci ^'^^^^ 
(3 = 0.005 meV [i^: the former excludes the ac-cycloidal 
order, while the latter generates an energy barrier be- 
tween ah- and 6c-cycloidal states. By minimizing the en- 
ergy, we obtained a stable structure of DW+a/+c- Then 
local electric polarization p was calculated on the basis 
of the inverse DM mechanism. As shown in Fig. 4, the 
calculated DW+a/-i-c structure is thick (~ 20 unit cells), 
reflecting the Heisenberg nature of the constituent spins. 



4 



Within the DW_|_q/_|.c, the cycloid plane continuously ro- 
tates from ab to be and thus the local polarization con- 
tinuously rotates from p||a to p||c. The latter feature is 
reminiscent of the ferromagnetic Neel wall in the Heisen- 
berg spin systems. 

The gradual rotation of polarization in the thick mul- 
tiferroic DW is a remarkable feature, in contrast to the 
situation in the conventional ferroelectric DWs, which 
is Ising-like and atomically thin 21, 22, [2^. Note that 



Rochelle salt, which show the weak polarization compa- 
rable to that of DyMnOa [2^ , is expected to exhibit the 
DW of a few unit cell thickness 2^, 2^- Therefore the 
thick multiferroic DW is attributable not merely to its 
small magnitude of polarization but to the ferroelectric- 
ity of magnetic origin. The motion of such Heiscnberg- 
like thick DWs is generally via magnon excitations with 
a small gap; thus, the magnon excitations and the 1/t 
of thick-DW motion are expected to diminish gradually 
toward low temperatures because of the small gap. In 
DyMnOa, the gradual slowing down of 1/tc obeying a 
power law [the inset of Fig. 2(d)] indicates that a gap, if 
any, is comparable with or even smaller than the present 
temperature range, 5-16 K: the gap may originate from 
the electromagnon-excitation gap at the P flop transi- 
tion. The possible small gap also implies thick DWs. 
Therefore the present experimental and numerical results 
are consistent, suggesting that the multiferroic D'W±a/±c 
is the Heisenberg-like thick DW rather than the Ising-like 
thin DW. This may explain why the present multiferroic 
DW is well mobile at low temperatures even though it 
works also as the ferroelectric DW. Finally, on the basis 
of the scaling behavior found by Catalan et al. [s^ , such 
thick DW implies large-size multiferroic domains as ob- 
serbed in the different kind of multiferroics BiFeOs [3l[. 

In summary, we have investigated the dielectric disper- 
sion of the GMC in DyMnOa over a wide frequency range 
and found that GMC is attributable not to the softened 
electromagnon but to the local motion of the multifer- 
roic DW between 6c-plane spin cycloid (P||c) and ab- 
plane spin cycloid (J'lla) domains. The relaxation rate 
of this DW motion holds high (~10^ s~^) even at 5 K. 
These results suggest that the multiferroic DW emerging 
at the polarization flop transition is the Hcisenbcrg-likc 
thick DW in contrast to the Ising-like thin DW in con- 
ventional fcrroelectrics. This difference may explain the 
mobile multiferroic DW even at low temperatures. 
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